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Melnikov Method for Autonomous Hamiltonians 

Clark Robinson 

Abstract. This paper presents the method of applying the Melnikov method 
to autonomous Hamiltonian systems in dimension four. Besides giving an 
application to Celestial Mechanics, it discusses the problem of convergence of 
the Melniltov function and the derivative of the Melnikov function. 



1. Introduction 

The original Melnikov method applies to periodically forced Hamiltonian sys- 
tems in dimension two (one degree of freedom). In this paper, our main empha- 
sis is to show how a modification of this method applies directly to autonomous 
Hamiltonian systems in dimension four (two degrees of freedom). These original 
unperturbed system needs to be completely integrable. By restricting to the energy 
surface, it is possible to convert the autonomous system in dimension four into a 
periodically forced Hamiltonian systems in dimension two and then apply the usual 
Melnikov method. It is our contention that this is not the natural method to treat 
this type of system. In fact the calculations are much more direct and natural 
when done directly with the autonomous system. Thus in this paper we advocate 
applying the Melnikov method directly to the autonomous Hamiltonian system. 
These ideas go back at least to the papers of P. Holmes and J. Marsden, [?] and 
[?] . Also see [?] and [?] . In Section 3 we discuss the main theorems in the theory of 
autonomous systems and their proofs. In Section 4 we indicate how these theorems 
can be applied to the circular restricted three body problem. 

Because we want to connect our treatment with the material for periodically 
forced systems, we start by describing the periodically forced situation using our 
approach in Section 2. 

This paper also points out the connection between the convergence of the im- 
proper integral giving the Melnikov function and the smoothness and dependence 
on parameters of the stable manifold of a periodic orbit. (See the discussion of 



Lemma 3.8 and the remarks which follow it.) This connection has also been made 
recently in the context of Celestial Mechanics by J. Casasayas, E. Fontich, and A. 
Nunes, [?!. 
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A final point made in this paper is the possibility of differentiating the Melnikov 
function by differentiating the integrand of its representation as a improper integral. 
See Theorem 3.5 in Section 3. 



2. Periodically force systems 

Consider a Hamiltonian : ^ M, with associated Hamiltonian vector field 
X = Xh- Next, consider the periodically force perturbation, 

X{z)+eY{z,T) 

T=l 

where r is considered as a periodic variable of some period T. We assume that X 
has a homoclinic orbit to a hyperbolic fixed point (in the z-space). Let 70 be the 
associated periodic orbit in the (z,r)-space, S, and 

Tq = {(zo,To) zq is homoclinic to 7} 

be the set of homoclinic orbits in the extended phase space S. 

Since 70 is hyperbolic, it can be continued for e 7^ to a periodic orbit 7^ in the 
extended phase space S. For a equal s and u, let W'^{'-fe,X + eY) be the stable and 
unstable manifolds of 7^ respectively. In order to determine the orbits on Tq to a 
particular orbit on W^lje^X + eY) for e ^ 0, we let A^(zo,ro) be a (one dimensional) 
transversal to Tq at a point (zq, tq) e Tq. Then for a = s,u, set 

r (zo, TO, e) = Wi^,, X + eY)n iV(,„,,„). 

In our development, we use strongly the fact that these stable and unstable 
manifold depend differentiably on the parameter and position on the manifold, so 
are smooth in the space 5, i.e., for a = s and u, C'^(zo,To,e) is a smooth function 
of all its variables. 

The Hamiltonian function is constant on Tq and has nonzero gradient, so it is 
a good measure of the displacement perpendicular to Tq. With this motivation, we 
define the Melnikov function by 

M(zo, To) = ^ [i? o C(zo, TO, e) - o C^(zo, TO, e)] l^^o- 

Thus, M is a measurement of the infinitesimal separation of the stable and unstable 
manifolds as measured by the Hamiltonian function H. By parameterizing the 
solutions of the equations and restricting to a transverse cross section, AI can be 
considered a function of a real variable (as it usually is). Notice that we use the 
fact that the stable and unstable manifolds depend smoothly on the parameter to 
define M. 

There are two theorems related to M. The first states that a nondegenerate 
zero of M corresponds to a transverse homoclinic orbit of 7^ for e 7^ but small 
enough. We r efer to [?], [?], or other standard references for this theorem. (Also see 



Theorem 3^ in the next section.) The second theorem relates M to an improper 
integral. We state this second theorem in order to emphasis the fact that with the 
proper statement, it can be calculated in any coordinate system (not just symplectic 
coordinates). 
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Theorem 2.1. 

/oo 

/oo 
-CXD 

We remark that the integral in equality (Q) is the usual one derived for the 
Melnikov function: the integrand is the scalar cross product (component in the 
k direction) or wedge product of the two vector fields. This integrand must be 
calculated in symplectic coordinates. The the integral in equality (Q) is valid in 
any set of coordinates, even those which are not symplectic. See [?] or [?] for a proof 
of the second form of the integral. Also in the discussion below of the autonomous 
case, we indicate a sketch of the proof why the above theorem is true. The proof in 
[?] uses the smoothness of the stable and unstable manifolds and not the uniform 
estimate for < t < oo on the difference of solutions for e 7^ and e = as many 
proofs do. (See e.g. [?] and [?].) 

3. Autonomous Hamiltonians in four dimensions 

In this section, the phase space is some four dimensional space which we label 
as S***: in Celestial Mechanics applications 5"* is usually or some product of 
copies of the circle (periodic variables) and copies of the reals. The autonomous 
Hamiltonian H'^ : S"' ^ R is given as 

i?^(z) =iJo(z)+ei?i(z) + 0(e2). 

The corresponding Hamiltonian vector field is given as 

X'-{z)^Xo{z)+eY{z)+0{e^). 

The basic assumptions are as follows. 
(Al) We assume that Xq is completely integrable, with a second independent 

integral F. (The function F is mostly likely not an integral of X'^ for e ^ 0.) 
(A2) We assume that Xo has a hyperbolic periodic orbit 70 with set of homoclinic 

orbits 

To -comp[l¥"(7o,^o)\7o] 
= comp[H/^(7o,Xo) \7o]. 

Here Fq is one of the components of the stable and unstable manifolds indi- 
cated and has dimension two because of the second integral. 

Example 3.1. The easiest example of the above type of system is one which 
decouples for e = into a subsystem with a homoclinic orbit and another That is 
oscillator: e.g., assume Jfo(z, w) = F{z) + G(w) where 

has a homoclinic orbit to the fixed point at 0, and 

G{wi,'W2) = ^{wl +wl) 
is an oscillator. A simple coupling is of the type Hi{z,w) = ZiWi. 
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Example 3.2. P. Holmes and J. Marsden [?] gave an example of a system 
which docs not decouple. (Compare this example with the closely related one from 
Celestial Mechanics given in [?].) Let 

2 7"2 

H'ip, q,e,I) = ^ + U{q) + ^ + esm{e), 

so 

i = -ecos{9) 

and F = I is a second integral for e = 0. The potential function U{q) must be 
chosen so that there is a homoclinic orbit in the {p, g')-subsystem for e = and / 
fixed. 

We set up the definition of the Melnikov function using transversals much as 
we did above for the periodically forced case. For each point zq G Tq, let N^o be a 
two dimensional plane perpendicular to Fg at Zq . Let 7^ be the continuation of the 
periodic orbit 70 for e 7^ with 7^ C {H'^)~^{ho) where ho = Ho(7o) is the energy 
of the unperturbed orbit. Then for a = s and u and e 7^ 0, let 

r(zo,/io,e) = A^zo nVF"(7e,X') 

be the continuation of Zq. Notice that both the stable and the unstable manifolds 
continue to lie in the same energy surface {H'^)~^{ho), so Hq is not a good measure 
of the separation; on the other hand, the second integral F is constant on Fq and 
is independent of Hq so it is a good measure of the distance perpendicular to Fq 
within this energy surface. The Melnikov function, M : Fq — > M is then defined as 
the infinitesimal displacement of the stable and unstable manifolds as measured by 
the second integral F: 

M(zo, ho) = §-^[Fo C"(zo, ho, e)-Fo C{zo, ho, e)] \^^^. 

By expressing M in terms of specific coordinates and restricting to a transverse 
section, it is possible to give it as a function of a single variable (plus the parameter 
ho)- See the example in the next section. 

Again, there are two theorems about the Melnikov function. The first states 
that a nondegenerate zero of M corresponds to a transverse homoclinic point. The 
second theorem indicates how to calculate M in terms of an improper integral. In 
general for this case, the integral representation of M is only conditionally conver- 
gent, so the times must be chosen correctly in the improper integral. The third 
theorem below justifies differentiating under the integral sign to calculate the de- 
rivative of the Melnikov function. 

Theorem 3.3. Make the basic assumptions A1-A2). If M{-,ho) : Fo R has 
a nondegenerate zero at Zq S Fq, then X'^ has a transverse homoclinic orbit z^ 
nearby for small enough nonzero e. 

Theorem 3.4. Make the basic assumptions (A1-A2). Assume thatTj and —T* 
are chosen so that the points (j)o{Tj, zq) and ^o{—Tj,zo) converge to the same point 
e 70- Then 

M(zo, ho) = lim / ' {DF ■ F)^„(t,zo) dt. 
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Theorem 3.5. Make the basic assumptions (A1-A2) as in Theorem \3.^ but 
assume that the integral converges absolutely (not conditional convergence). (See 
remarks following Lemma \3.^ for conditions which imply this assumption.) Then 



/OO 
D,{{DF-Y)^^^^t,^))vdt. 
-oo 



In the above equation, Dz is differentiation with respect to the initial spatial coor- 
dinate z, and is calculated in the direction of the vector v. 



Remark 3.6. The integral in Theorems 3.4 and 3.5 can be calculated in any 
coordinates, even non-symplectic coordinates. This is relevant to applications in 
Celestial Mechanics because McGehee coordinates are non-symplectic coordinates. 
See Section 4. 

Remark 3.7. Although the convergence is only conditional in general, there 
are several circumstances in which the convergence is absolute. We make this more 
explicit below. 



The proof of Theorem 3.3 follows as usual from the Implicit Function Theorem, 
so we do not give any details. See [?]. 



Idea of the proof of Theorem |3^. As in the usual proof, 
M(zo,/io)- {DF ■Y)^^^(t^^^)dt + R{-T;,T,) 

J -T' 

3 

where the remainder term 

Ri-T;,T,) =^ [F o C(zo, e), e)] 

-^[^o0(r„C(zo,e),6)]|^^„. 

In [?], we prove a lemma about the convergence of the remainder. The statement 
of Lemma 2.1 in [?] only includes the final statement of the following lemma but 
the proof in that paper implies the following extended statement. Clearly Lemma 



3.8 implies the Theorem 3.4, □ 



Lemma 3.8. Assume thatTj and ~T* are chosen so that c/jq^—T* ,z,o) converges 
to z^ and (t>o(Tj,zo) converges to z^ as j ^ oo. Then 

^.[Fo0(-i;*,C(zo,e),e)]L^o 



converges to 



and 



converges to 



de 

d_ 

We 



^[^o^(T„C(zo,6),e)]|^^„ 



l[Fo^{zt,e)]l^^. 
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Here 7(2, e) is the Junction which gives the perturbed orbit 7^ Jar e ^ as a graph 
over the unperturbed orbit. In particular, if = then R(—T*,Tj) goes to zero 
as i goes to infinity. 



Remark 3.9. By Lemma |3.8| , the improper integral in Theorem 3.4 always 
gives M{zq, Hq) provided Tj and —T* are chosen so that = z+ . 

Remark 3.10. If DF(^z) = for all z G 70, then the convergence is not condi- 
tional and the integral is just the improper integral from —00 to 00. This feature 
is essentially what is true in the proof for the usual Melnikov function when using 
the function H which vanishes at the fixed point. 

Remark 3.11. If the periodic orbit does not move with e, = 7o, then the 
convergence is also a usual improper integral (and not conditional) since the limiting 
values are zero and do not need to cancel each other. When applying this theorem 
to parabolic orbits which are asymptotic to periodic orbits at infinity, this is often 
the case since the orbits at infinity are not affected by the small coupling terms. 
See the example of the circular restricted three body problem in Section 4. 



Remark 3.12. What makes the convergence work in the proof of Lemma 3. J 
is the stable manifold theory. Casasayas, Fontich, and Nunes proved in [?] that this 
type of result is often true even when the periodic orbit is weakly hyperbolic (from 
higher order terms). Notice that they emphasize that it is possible to prove a stable 
manifold theorem in the parabolic case with a particular smoothness. Then then 
prove that this smoothness of the stable manifold is sufficient to imply a convergence 
of the above type to give the Melnikov function as an improper integral. They treat 
the usual time periodic perturbations but their results go over to the above context 
of autonomous systems as well. 



Proof of Theorem p.5| . The fact that we can differentiate under the inte- 
gral sign follows from a theorem in analysis about differentiating improper integrals. 
What is needed to differentiate under the integral sign of an improper integral is 
that the integrals for M(z, /iq) and -Dz-A^(z./io)'^ converge uniformly in z. See [?] 
or [?] . In our situation the convergence of both integrals is uniform for the same 



reason that the integral converges in the proof of Lemma 3.8 as given in [?]. □ 



4. Planar Restricted Three Body Problem Example 

Both J. Llibre and C. Simo in [?] and Z. Xia in [?] applied the usual Melnikov 
method to the circular planar three body problem. In this section, we indicate 
how the same results can be obtain (possibly) more directly by using the method 
for the autonomous systems indicated in the last section. Our remarks about the 
applicability in non-symplectic coordinates also justifies making the calculations in 
McGehee coordinates. 

Below we sketch the introduction of the coordinates into the problem. See 
the above references for more detailed explanations and treatment. Our treatment 
most closely follows that in [?]. 

In this situation, the third mass TO3 = 0, mi = is a small parameter, and 
TO2 = 1 — /i. Jacobi coordinates are used as a starting place. The problem is circular 
so the variable 

Q = (cos(i), sin(t)) 
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gives the motion of the vector from m2 to mi . The position of ms relative to mi 
and m2 is given by q. 

Next, McGehee coordinates are introduced. The variable x measures the dis- 
tance from infinity, 

= |q|, 

SO X goes to zero as |q| goes to infinity. The direction to ma is given by 6, 

7-7 = (cos^, sin^). 

|q| 

This angle is measured in rotating coordinates by introducing the 27r-periodic vari- 
able s — t — 6. The momentum variable for m^, p is broken down into radial 
and angular components by introducing the variables y and p respectively by the 
equation 

p = [y + x^pi\ [cos ^ + « sin ^] . 

(This last equations tiscs the notation of complex variables for points in M^.) The 
energy equation in these variables becomes 

H>^ = \y^ + \x'p^-U>^ 

where = . 

When = 0, p = 0, so p, which is the angular momentum of TO3 (the Jacobi 
constant), is a second integral. After a change of the time variable, the equations 
become 

1 3 

X = -^x y 

y' = -x'^ + x^p^ + pgi{x, s) + 0{n'^) 
p' = pg2{x,s) +0{p?) 
s' = 1- x^p. 

Here s is considered a 27r periodic variable. 

When p, = 0, = ^jj^ is a fixed point in {^y^ -space. This fixed point 

corresponds to a periodic orbit in the total space (with s the periodic variable). 
This periodic orbit is weakly hyperbolic but it has smooth stable and unstable 
manifolds which correspond to the parabolic orbits of the original problem. We 
parameterize the homoclinic orbits 

X = ^o{t,po,so) 

y = v{t,po,so) 

s = so + cr{t, po) 

where we take the initial conditions so that 77(0, pa, sq) = 0, i.e. the orbits cross the 
transverse section y = at t = 0. Fixing po at t = is like fixing the energy of the 
system. 
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For fi = 0, p is constant and so is a second integral of the system. Thus the 
Melnikov function is 



M (so, Po) = / (Dp) ■ i^0o(t,Po,so) 

92 o(j)o{t,po,sa)dt 
Co(i) sin(so + a{t))S{t, po, sq, ) dt, 



where 

The above form of the Melnikov function is derived in [?] by starting with the form 
arising from a periodically forced system and then integrating by parts. This type 
of derivation is not surprising since the equivalence of tbe usual Melnikov integral 
and that obtained for autonomous sytems is proved in [?] using integration by parts. 
The point that is being made here is the above treatment is very direct and seems 
less ad hoc then what is found in [?]. (I am certain that Z. Xia was very aware of 
the representation of the Melnikov integral by means of a second integral and that 
he just preferred using the more familiar form.) Notice that the orbit at infinity 
does not vary with p so the integral converges absolutely. 

Once the above form is obtained, symmetry considerations imply that there is 
a zero of M for so = and tt, M(0, po) = M(7r, po) = 0. 

The orbit at infinity does not vary with p, so the above integral representation 

dM 



of M is not conditional. Therefore by Theorem 3^, — — (so,/0o) can be calculated 

aso 



by differentiating under the integral sign with respect to so which appears in a 
very simple way. (The integral of this derivative with respect to so also converges 
uniformly in sq.) Special arguments are needed to get that the improper integral 
dM 

representing — — (so;Po) for sq = or tt is nonzero, and so the corresponding zero 

of Af(-,po) is a nondegenerate zero as a function of so. These arguments which 
verify this fact done very differently in [?] and [?] . 
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